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ABSTRACT 

This  report  presents  an  evaluation  of  tire  steady-state  capability  of  the  turbulence  models 
available  in  the  commercial  CFD  code  FLUENT  6.0,  for  their  application  to  the  simulation  of 
hydrofoil  turbulent  boundary  layer  separation  at  high-Reynolds  numbers.  Four  widely 
applied  two-equation  RANS  turbulence  models  were  qualitatively  and  quantitatively 
assessed  through  comparison  with  high-quality  experimental  data  at  Reynolds  numbers  of 
8.284 '106  and  1.657'107.  The  turbulence  models  evaluated  were  the  Standard  k-e  model,  the 
Realisable  k-e  model,  the  Standard  k-w  model  and  the  Shear-Stress-Transport  (SST)  k-w 
model.  It  was  concluded  that  the  Realisable  k-e  turbulence  model  used  with  enhanced  wall 
functions  and  near-wall  modelling  techniques,  consistently  provided  superior  performance  in 
predicting  the  hydrofoil's  performance  parameters  and  flow  characteristics. 
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Hydrodynamic  Flow  Simulations 


Executive  Summary 

Hydrodynamic  flows  of  interest  to  designers  of  large  surface  and  sub-surface  marine 
vessels  and  their  propulsors  are  typically  characterised  by  high-Reynolds  numbers  and 
usually  some  form  of  flow  unsteadiness.  To  numerically  model  such  flows  poses 
significant  challenges  to  the  current  range  of  'two-equation'  RANS  (Reynolds- 
Averaged  Navier-Stokes)  turbulence  models. 

In  an  effort  to  determine  the  most  accurate  turbulence  model  for  high-Reynolds 
number  hydrodynamic  flows,  an  evaluation  of  the  steady-state  capability  of  the 
turbulence  models  available  in  the  commercial  CFD  code  FLUENT  6.0  has  been 
undertaken.  This  evaluation  is  limited  to  applications  to  the  simulation  of  hydrofoil 
turbulent  boundary  layer  separation  at  high-Reynolds  numbers. 

Four  widely  applied  two-equation  RANS  turbulence  models  were  qualitatively  and 
quantitatively  assessed  through  comparison  with  high-quality  experimental  data  at 
Reynolds  numbers  of  8.284x106  and  1.657xl07.  In  this  preliminary  investigation,  the 
turbulence  models  evaluated  were  the  Standard  k-s  model,  the  Realisable  k-s  model, 
the  Standard  k-co  model  and  the  relatively  new  Shear-Stress-Transport  (SST)  k-co  model. 

It  was  concluded  that  the  Realisable  k-s  turbulence  model  used  in  conjunction  with 
enhanced  wall  functions  and  near-wall  modelling  techniques,  consistently  provided 
superior  performance  in  predicting  the  hydrofoil's  performance  parameters  and  flow 
characteristics.  This  turbulence  model  and  wall  treatment  accurately  predicted  the 
pressure  coefficient  distribution,  the  lift  and  drag  coefficients,  the  pressure  surface 
boundary  layer  characteristics,  the  separation  of  the  suction  surface  boundary  layer 
and  the  flow  features  in  the  near  wake  region. 
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1.  Introduction 

Hydrodynamic  flows  of  interest  to  designers  of  large  surface  and  sub-surface  marine 
vessels  and  their  propulsors  are  typically  very  complex.  The  length  scales  and  the 
relative  flow  velocities  in  both  the  free-stream  and  induced  fields  of  these  vessels 
produce  flows  characterised  by  very  high  Reynolds  numbers.  Together  they  pose  the 
greatest  challenge  to  both  the  experimental  and  theoretic  researches  in  modem  fluid 
dynamics.  Such  flows  can  typically  exceed  Reynolds  numbers  of  107  [1],  yet  can  reach 
magnitudes  of  109  and  1010  [2].  Secondly,  these  flows  exhibit  of  high  unsteadiness  or 
time-dependent  phenomena,  such  as  the  shedding  of  vortices  from  propeller  blades, 
control  surfaces  and  other  submerged  appendage.  Thirdly,  with  the  occurrence  of 
cavitation  these  flows  can  also  be  multiphase. 

The  expectations  placed  upon  solution  methodology,  computational  software  and 
computer  hardware  to  simulate  these  hydrodynamic  flows  using  numerical  methods 
are  immense.  With  the  emergence  of  parallel  computing  and  its  provision  of  reduced 
global  memory  and  run-times  through  multiple  processors,  numerical  solutions  are 
becoming  more  achievable  for  increasingly  complex  problems  [3].  Parts  of  larger  three- 
dimensional  flow  problem  that  were  previously  simulated  individually  are  now  being 
simulated  as  complete  configurations,  such  as  full  submarine  geometries  comprising 
hull,  sail,  stern  appendage,  propeller  and  moving  control  surfaces  for  both  steady  and 
unsteady  flow  scenarios  [3]. 

Given  that  both  Direct  Numerical  Simulation  (DNS)  and  Large-Eddy  Simulation  (LES) 
are  still  untenable  approaches  in  practical  flows,  the  modelling  of  most  practically 
relevant  turbulent  flows,  such  as  most  hydrodynamic  flows,  is  widely  based  upon  the 
solution  of  Reynolds- Averaged  NS  (RANS)  equations.  As  a  result,  accurate  modelling 
of  Reynolds  stresses  becomes  crucial.  There  are  various  turbulence  modelling 
approaches.  Questions  remain  about  whether  or  not  the  available  turbulence  models 
predict  the  correct  flow  physics  and  which  turbulence  models  are  best  suited  to  what 
type  of  flows.  It  is  agreed  [3]  that  it  is  both  naive  and  dangerous  to  assume  the  validity 
of  a  turbulence  model  over  a  wide  range  of  flows.  It  is  also  agreed  that  unique 
modification  or  refinement  of  existing  turbulence  models  for  specific  flow  applications 
may  provide  improved  performance.  Many  evaluations  of  two-equation  RANS 
turbulence  models  documented  in  the  literature  have  been  undertaken  in  order  to 
identify  the  most  suitable  model  for  the  particular  flow  of  interest. 

Speziale  et.  al.  [4]  applied  a  variety  of  two-equation  (k-s  and  k-co  variants)  models  and 
their  k-r  model  to  near- wall  turbulent  flows.  They  concluded  that  due  to  the 
excessively  high  levels  of  destruction-dissipation  and  turbulence  viscosity  near  the 
wall,  the  predictions  of  the  k-s  models  tested  were  severely  compromised.  As  far  as  the 
k-co  model  of  Wilcox  et.  al.  [5,  6]  was  concerned,  this  model  had  the  same  constraints  as 
the  k-s  models,  yet  returned  excessive  turbulent  kinetic  energy  near  the  wall.  However, 
for  many  applications,  the  k-co  model  provided  adequate  prediction  of  the  skin  friction 


1 


DSTO-TR-1564 


at  the  wall  and  the  turbulence  statistics  away  from  the  wall.  Based  on  the  numerical 
tests  of  Speziale  et.  al.  [4],  their  k-t  turbulence  model  seemed  to  have  improved 
predictions  within  the  turbulent  boundary  layer. 

Lien  &  Leshziner  [7]  examined  both  a  non-linear  eddy-viscosity  model  and  a  second- 
moment  closure  model  by  applying  them  to  2D  separated  flow  from  a  high-lift  airfoil. 
They  concluded  that  the  second-moment  closure  produced  superior  predictions 
relative  to  both  the  non-linear  eddy-viscosity  model  and  the  realisable  Reynolds  stress 
algebraic  model  of  Shih  et.  al.  [8].  The  second-moment  closure  model  had  the  inherent 
ability  of  representing  the  interaction  between  streamline  curvature  and  turbulence 
anisotropy,  and  the  effect  of  anisotropy  on  shear  stress.  Nevertheless,  the  second- 
moment  model  failed  to  predict  the  onset  of  separation  at  the  correct  angle  of  attack. 
The  failure  was  attributed  to  the  excessively  high  levels  of  shear  stress  near  the  wall, 
preventing  separation. 

Pajayakrit  &  Kind  [9]  evaluated  four  turbulence  models  for  predicting  the  development 
of  plane  and  curved  wall  jet  flows.  The  models  evaluated  were  the  Baldwin  &  Lomax 
[10]  algebraic  model,  the  Dash  et.  al.  [11]  k-e  model,  the  Wilcox  [5]  k-co  and  the  Wilcox 
[6]  multi-scale  model.  Without  any  modification,  all  models  demonstrated  significant 
deficiencies,  although  the  Dash  et.  al.  [9]  k-e  model  performed  reasonably  well  on  an 
overall  sense,  particularly  for  the  curved  wall  jet  flow.  Fine  tunning  of  the  Wilcox  [5,6] 
models  was  undertaken  by  adjusting  the  model  constants,  however,  improvement  was 
achieved  by  the  modified  k-co  model  only. 

Unlike  the  work  presented  by  Speziale  et.  al.  [4],  Lien  &  Leshziner  [7]  and  Pajayakrit  & 
Kind  [9],  Kim  et.  al.  [12]  evaluated  three  widely  applied  and  economic  turbulence 
models  using  a  practical  problem;  the  hull  form  design  of  a  very  large  crude-oil  carrier. 
The  turbulence  models  evaluated  were  the  Standard  k-e  model,  the  RNG  k-e  model  and 
the  Realisable  k-e  model.  Emphasis  was  placed  on  their  ability  to  accurately  predict  the 
nominal  wake  distribution  at  the  propeller  plane,  a  region  of  strong  bilge  vortices  and 
complex  flow  behaviour,  as  flow  information  within  this  region  is  vital  for  propeller 
design.  It  was  determined  that  the  Realisable  k-e  turbulence  model  accurately  predicted 
the  strength  and  location  of  the  bilge  vortex,  however,  both  the  Standard  k-e  and  RNG 
k-e  turbulence  models  failed  to  predict  this  flow  feature. 

This  report  is  not  intended  as  an  authoritative  reference  for  turbulence  model 
applications  for  all  hydrodynamic  flow  simulations,  nor  have  we  aimed  to  discuss  the 
strengths  and  weaknesses  of  the  turbulence  models  at  a  detailed  level.  Rather,  like  the 
work  presented  by  Kim  et.  al.  [12],  it  report  is  intended  as  a  guide  to  researchers  and 
designers  who  are  working  on  a  similar  hydrodynamic  problem. 

The  hydrodynamic  problem  of  interest  in  this  report  is  the  numerical  simulation  of  the 
turbulent  hydrofoil  flow  at  high  Reynolds  numbers,  where  the  boundary  layer  on  the 
suction  surface  is  subjected  to  a  large  adverse  pressure  gradient.  This  category  of 
turbulent  flow  poses  major  modelling  challenges  for  the  predictive  capabilities  of  any 
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turbulence  model.  Such  turbulent  boundary  layer  separation  occurs  at  the  trailing  edge 
of  hydrofoils  operating  at  high  Reynolds  numbers,  resulting  in  the  formation  of 
unsteady  oscillating  wakes.  Under  certain  conditions  the  vortex-shedding  period  of  the 
wake  may  become  coincident  with  the  resonant  frequency  of  the  hydrofoil  or  parts  of 
it,  resulting  in  the  emission  of  a  high-level  tonal  noise.  Such  an  emission  from  lifting 
surfaces  associated  with  naval  vessels  is  highly  undesirable,  as  this  noise  is  a  dominant 
source  of  information  for  underwater  detection  systems.  This  preliminary  study 
evaluates  the  steady-state  capability  of  several  two-equation  RANS  turbulence  models 
using  the  commercial  CFD  code  FLUENT  6.0.  The  four  turbulence  models  evaluated 
are  the  Standard  k-s,  the  Realisable  k-s,  the  Standard  k-co  and  the  relatively  new  Shear- 
Stress-Transport  (SST)  k-co,  as  detailed  in  the  literature  [13]. 

Also,  in  an  effort  to  provide  researchers  with  a  guide  to  the  application  of  wall 
functions  in  such  a  flow  scenario,  three  different  wall  treatments  are  investigated  with 
respect  to  their  suitability.  These  include  two  different  wall  functions,  the  standard  and 
non-equilibrium  variants,  and  enhanced  wall  treatment  employing  near-wall 
modelling  techniques. 

The  developed  CFD  model  is  based  on  the  characteristics  and  parameters  of  the  recent 
experimental  research  into  high-Reynolds  number  hydrofoil  flow,  undertaken  by 
Bourgoyne  et.  al.  [1,14  and  15],  allowing  direct  comparison  to  high-quality 
experimental  data. 
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2.  Problem  Overview 


The  main  objective  of  the  experimental  research  conducted  by  Bourgoyne  et.  al.  [1, 14 
and  15]  was  to  identify  and  document  the  dominant  flow  features  in  the  trailing  edge 
region  of  a  hydrofoil  operating  at  chord-based  Reynolds  numbers  approaching  10s.  In 
addition  this  research  has  provided  a  database  of  high-quality  experimental  data  for 
the  development,  refinement  and  validation  of  turbulence  models  for  high-Reynolds 
number  applications. 

The  research  was  conducted  in  the  William  B.  Morgan  Large  Cavitation  Channel  (LCC) 
in  Memphis,  USA.  The  test  hydrofoil  was  fully  two-dimensional,  spanning  the  test 
section  (3.05m),  with  a  chord  length  (C)  of  2.134m.  The  cross-section  profile  was 
representative  of  a  generic  naval  propeller  of  moderate  thickness  ( t )  and  camber  (/), 
utilising  a  NACA-16  airfoil  profile  (t/C  =  0.08,  f/C  =  0.032)  with  two  modifications. 

The  first  modification  was  that  the  pressure  surface  of  the  hydrofoil  was  flat  beyond 
28%  chord  length  continuing  to  the  trailing  edge.  The  second  modification  was  that 
near  97%  chord  length,  the  suction  surface  steeply  terminated  at  the  trailing  edge  at  an 
angle  of  44°  to  the  flat  pressure  surface,  a  characteristic  typical  of  anti-singing  trailing 
edges.  A  detailed  diagram  of  the  hydrofoil  geometry  is  illustrated  in  Figure  1  and  the 
anti-singing  trailing  edge  geometry  is  detailed  in  the  literature  [1], 

Experiments  were  conducted  at  flow  velocities  of  0.5m/s  to  18.3m/s,  corresponding  to 
chord-based  Reynolds  numbers  of  1.4x106  to  6xl07.  The  hydrofoil  was  aligned  with  the 
flow  (i.e.  zero  angle  of  attack)  and  its  mounting  method  within  the  water  tunnel 
produced  good  span-wise  two-dimensional  flow  [1],  Laser-Doppler-Velocimetry  (LDV) 
and  Particle-Imaging-Velocimetry  (PIV)  measurement  techniques  were  used  to 
determine  mean  flow  velocities  and  turbulence  statistics  from  the  separating  boundary 
layer  flow  near  the  hydrofoil's  trailing  edge.  Static  pressure  measurements  on  the 
surface  of  the  hydrofoil  were  made  using  seventeen  suction-side  and  thirteen  pressure- 
side  taps. 

The  mean  flow  velocities  at  2500  co-ordinate  positions  within  the  trailing  edge  region 
were  measured  using  the  LDV  system.  A  nominal  sampling  period  of  36  seconds  per 
co-ordinate  was  used,  with  the  number  of  samples  varying  between  500  and  12000.  The 
vortex  shedding  period  of  the  hydrofoil  was  theoretically  calculated  to  be  60  ms  at  a 
flow  speed  of  3  m/ s  and  10  ms  at  18.3m/s.  Therefore,  the  sampling  period  was  equal  to 
between  600  to  3600  shedding  periods,  allowing  accurate  resolution  of  die  mean 
properties  within  the  unsteady  flow.  The  LDV  system  was  determined  to  have  a 
maximum  calibration  bias  range  in  the  normalised  mean  velocities  of  ±0.02  at  3m/s  to 
±0.006  at  18.3m/s  [1].  The  positioning  error  of  the  LDV  system  was  assessed  as 
±0.1mm,  equivalent  to  ±0.0047%  of  the  chord  [1]. 
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Two  main  conclusions  were  drawn  from  the  experimental  testing.  Firstly,  the 
hydrofoil's  near  wake  flow  features  appeared  to  be  Reynolds  number  dependent  and, 
secondly,  the  observation  of  Reynolds  number  dependency  is  consistent  with  suction- 
side  boundary  layer  separation  occurring  closer  to  the  trailing  edge  with  increasing 
Reynolds  number. 

Although  the  application  of  numerical  methods  to  hydrodynamic  lifting  surface  flow  is 
widespread  in  the  literature  [16, 17, 18  &  19]  and  has  greatly  matured,  there  are  specific 
challenges  posed  to  numerically  modelling  the  experimental  testing  conducted  by 
Bourgoyne  et.  al.  [1,  14  &  15].  Firstly,  the  experimental  testing  was  conducted  at  a 
Reynolds  numbers  well  above  the  range  at  which  the  turbulence  models  were 
validated.  Secondly,  the  numerical  model  must  accurately  predict  tire  separation  of  the 
turbulent  boundary  layers  from  both  the  pressure  and  suction  surfaces  at  the  trailing 
edge.  Separation  itself  is  a  highly  complex  physical  phenomenon,  and  it  also  introduces 
mesh  issues  in  numerical  modelling.  Thirdly  and  most  importantly  for  hydroacoustic 
research,  the  numerical  model  must  also  be  able  to  correctly  resolve  the  shedding  of 
vortices  and  the  formation  of  an  oscillating  wake. 

Of  particular  importance  to  this  preliminary  steady-state  analysis  is  boundary  layer 
separation.  In  this  flow  configuration  with  a  slender  and  streamlined  body,  the 
numerical  modelling  of  turbulent  boundary  layer  separation  is  more  challenging  due 
to  its  Reynolds  number  dependency.  It  has  been  shown  experimentally  by  both 
Baubeau  &  Latorre  [17]  and  Bourgoyne  et.  al.  [1]  that  at  a  constant  angle  of  attack,  a 
turbulent  boundary  layer  resists  separation  further  along  the  surface  of  a  hydrofoil, 
separating  closer  to  the  trailing  edge  with  increasing  Reynolds  number. 

Correct  modelling  of  turbulent  boundary  layer  separation  is  crucial  for  the  accurate 
prediction  of  the  coefficients  of  lift  and  drag.  Incorrect  prediction  of  the  separation 
location  can  lead  to  errors  in  the  magnitude  of  the  low-pressure  recirculation  bubble  at 
the  trailing  edge,  hence  over-  or  under-prediction  of  the  pressure  drag  of  the  hydrofoil. 

As  in  this  application  the  pressure  surface  of  the  hydrofoil  is  flat,  separation  of  the 
pressure  surface  boundary  layer  at  the  trailing  edge  tip  is  inevitable.  Therefore 
emphasis  is  placed  upon  the  prediction  of  the  suction  surface  boundary  layer 
separation  from  the  anti-singing  trailing  edge. 
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3.  Theory 

3.1  Governing  Equations 

The  governing  equations  for  the  turbulent  incompressible  flow  encountered  in  this 
research  are  the  steady-state  Reynolds-averaged  Navier-Stokes  (RANS)  equations  for 
the  conservation  of  mass  and  momentum.  In  the  literature  [13]  they  are  presented  in 
the  following  forms: 

Continuity: 

|-(w)=0  (!) 

Momentum: 


d 


dx. 


dp  d 

dx,  dXj 


(2) 


In  Equations  (1)  and  (2),  p  is  mean  density,  p  is  mean  pressure,  p  is  tire  molecular 
viscosity  and  -  pup  j  is  the  Reynolds  stresses.  To  correctly  account  for  turbulence, 

Reynolds  stresses  are  modelled  in  order  to  achieve  closure  of  Equation  (2).  The  method 
of  modelling  employed  utilises  the  Boussinesq  hypothesis  to  relate  the  Reynolds 
stresses  to  the  mean  velocity  gradients  within  the  flow.  Therefore  the  Reynolds  stresses 
are  given  by: 


pu,u  =  p, 


^  du  t  diij 
,  dXj  dx, 


f  ,  diii  ^ 

pk  +  p  — 

dx,  j 


(3) 


In  Equation  (3),  pt  is  the  turbulent  (or  eddy)  viscosity  and  k  is  the  turbulent  kinetic 
energy.  For  two-equation  turbulence  models  such  as  the  k-s  and  k-co  variants,  the 
turbulent  viscosity  is  computed  through  the  solution  of  two  additional  transport 
equations  for  the  turbulent  kinetic  energy,  and  either  the  turbulence  dissipation  rate,  e, 
or  the  specific  dissipation  rate,  to. 
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3.2  Two-Equation  Turbulence  Models 

3.2.1  Standard  k-s 

The  Standard  k-e  turbulence  model  is  presently  the  most  widely  applied  turbulence 
model  to  practical  engineering  flows  as  it  is  robust,  economical  and  provides 
reasonable  accuracy  for  a  wide  range  of  flows.  The  transport  equation  for  k  is 
physically  correct,  however  the  transport  equation  for  e  is  heavily  modelled  [13].  The 
modelled  transport  equations  for  k  and  s,  for  steady-state  and  incompressible  flow,  are 
given  in  the  literature  [13]  as: 


In  Equations  (4)  and  (5),  ok  and  ae  are  the  turbulent  Prandtl  numbers  for  k  and  s 
respectively.  The  production  of  turbulence  kinetic  energy,  Gk,  is  approximated  in  a 
manner  consistent  with  the  Boussinesq  hypothesis  by: 

G,=fi,S2  (6) 

S  is  the  modulus  of  the  mean  rate-of-strain  tensor,  defined  by: 

S=^s~ 

„  if  duj  6ui 

S  = - -  + - 

J  2  ^  dxj  dxj  ^ 

The  dissipation  of  this  turbulence  kinetic  energy,  Yk,  is  defined  by: 

Yu  =  P£  (7) 

The  production  and  dissipation  of  s,  G,  and  Ye  respectively,  are  defined  as: 

C«=C„f(G1)  (8) 

k 
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Y,:  =  Ch:PT  (9) 

In  Equations  (8)  and  (9),  Ci^  and  €2,-.  are  constants.  The  turbulent  viscosity  is  computed 
through  a  formulation  of  k  and  e,  given  as: 


V, 


(10) 


In  Equation  (10),  Cf,  is  a  constant,  and  this  constant  and  the  remaining  constants 
applied  in  the  Standard  k-e  turbulence  model  are  equal  to: 

Cu  =  1.44,  C2t  =  1.92,  C„  =  0.09,  o*  =  1.0,  as  =  1.3 


3.2.2  Realisable  k-e 

The  Realisable  k-e  turbulence  model  differs  from  the  Standard  k-e  model  in  two 
important  ways.  Firstly  it  contains  a  new  formulation  for  the  turbulent  viscosity,  and 
secondly,  a  new  transport  equation  for  e  has  been  derived  from  an  exact  equation  for 
the  transport  of  the  mean-square  vorticity  function  [13] . 

In  the  Realisable  k-e  turbulence  model,  the  transport  equation  for  k  is  identical  to  that  of 
the  Standard  k-e  model  (Equation  (4)).  However,  the  production  and  diffusion  terms  in 
the  transport  equation  for  e  differ  slightly  and  are  defined  as: 

Gr  =  pCxSe  (11) 

y^pc2T±~T  (12) 

k  +  y/ve 

In  Equations  (11)  and  (12),  C2  is  a  constant,  while  Cj  is  defined  as: 

Cj  =  max  0.43,—— 

L  r/  +  5_ 


£ 

Similarly  to  the  Standard  k-e  turbulence  model,  the  turbulent  viscosity  is  computed 
using  Equation  (10),  however,  Cf,  is  no  longer  a  constant.  Cfl  is  calculated  using  the 
following  equations: 
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c„=- 


A  +  4 


kU* 


(13) 


Where: 


ss=i 


^  dui 
KdxJ 


duj 
dx . 


In  Equation  (13),  Ao  is  a  constant  and  the  remaining  variable.  As,  is  calculated  using  the 
following: 


As  =  y[6  COS  (j) 


Where: 


SjjSjkSk, 


The  constants  applied  in  the  Realisable  k-s  turbulence  model  are  equal  to: 


ok  =  1.0,  =  1.2,  C2  =1.9,  Ao  =4.04 


3.2.3  Standard  k-co 

The  Standard  k-co  turbulence  model  incorporates  modifications  for  low-Reynolds 
number  effects  and  shear  flow  spreading.  The  shear  flow  spreading  rates  are  in  close 
agreement  with  measurements  for  far  wakes,  mixing  layers,  and  plane,  round,  and 
radial  jets,  and  is  thus  applicable  to  wall-bounded  and  free  shear  flows  [13]. 

The  modelled  transport  equation  for  k  in  the  Standard  k-co  turbulence  model  is  very 
similar  to  that  of  the  k-£  based  models.  The  production  of  turbulence  kinetic  energy  is 
defined  identically  to  that  of  the  k-s  models  (Equation  (6)),  however  the  Standard  k-co 
model  has  a  different  formulation  for  the  diffusion  of  k,  defined  as: 

r*  =  Pphfp  koi  (14) 


Where  /T  is  a  constant  and: 
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fr=\  1  +  680  xl 

1  +  400 x; 


Xk^  0 

Xk  >0 


_  1  dk  dco 
co 3  dxj  dXj 

With  replacement  of  e  for  co  in  Equation  (5),  the  transport  equation  for  co  in  the 
Standard  k-co  model  is  created  from  that  of  a  in  the  k-e  based  models.  However,  the 
production  and  diffusion  terms  of  co  are  significantly  different  to  those  of  e.  The 
production  and  diffusion  of  co  in  the  high-Reynolds  number  form  of  the  Standard  k-co 
model  are  defined  respectively  as: 


G-IG‘ 


Wire  re  /?  is  a  constant  and: 


Y,„  =  pPifpV2 


l  +  70Zft 


l  +  80jft 


t b> y 

Again,  another  deviation  of  the  Standard  k-co  from  the  previous  turbulence  models  is 
the  computation  of  die  turbulent  viscosity.  For  die  high-Reynolds  number  form  of  the 
Standard  k-co  model,  the  turbulent  viscosity  is  calculated  using: 

f,=—  (17) 

CO 

The  constants  applied  in  the  high-Reynolds  number  form  of  the  Standard  k-co 
turbulence  model  are  equal  to: 


ok  =  2.0,  cr,„  =  2.0,  pi  =  0.09,  /?,  =  0.072 
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3.2.4  Shear-Stress-Transport  (SST)  k-co 

The  SST  k-co  turbulence  model  is  a  conglomeration  of  the  robust  and  accurate 
formulation  of  the  Standard  k-co  model  in  the  near-wall  region,  with  the  Standard  k-£  in 
the  far  field.  The  SST  k-co  is  more  accurate  and  reliable  for  a  wider  class  of  flows  than 
the  Standard  k-co,  including  adverse  pressure  gradient  flows  [13]. 

In  the  SST  k-co  turbulence  model,  the  modelled  transport  equation  for  k  is  identical  to 
that  of  the  Standard  k-co  model,  however,  the  transport  equation  for  ca  differs  slightly 
with  the  inclusion  of  a  cross-diffusion  term,  D,„. 


The  production  of  k  in  the  SST  k-co  model  is  defined  in  the  same  manner  as  per  the 
Standard  k-co  model,  hence  the  k-e  type  models,  by  Equation  (6).  However,  variations 
exist  in  the  production  of  co  with  comparison  to  the  Standard  k-co  model.  The 
production  of  co  is  given  by: 

Gm^Gk  (18) 

K 

Where: 


a 


00,1 


a„,2  = ' 


Pu 2 

p: 


The  dissipation  of  the  k  and  co  are  defined  differently  to  the  Standard  k-co  model  by  the 
following  respectively: 


Where: 


Yk  =  pfi'kco 

A  =  *iAi+(i  -*iK2 


(19) 

(20) 


(21) 
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In  Equation  (21),  /?,,  and  /?,  2  are  constants  while  the  Fi  is  determined  from  the 
following: 

Fx  =  tanh(<&?) 


O,  =  min  max 


4k  500/yl  4  pk 


0.09  ay’  0’2o)  j  G^Dy- 


n+  -11  dk  dco  —  _2o 

D„=  max  2 p - — —  ,10 

crlo  2  co  CXj  dxj 


In  the  previous  equations,  y  is  the  distance  to  the  next  surface  and  D*o  is  the  positive 
portion  of  the  cross-diffusion  term  (to  be  introduced  at  a  later  stage). 

Unlike  the  Standard  k-co  model  where  the  turbulent  Prandtl  numbers  for  k  and  co  (<n 
and  cr,.,)  are  constants,  in  the  SST  k-co  model  they  are  non-constants  defined  by: 


F?  is  defined  using  the  following: 


rA\l  <Ti t.2 


Ex  ,  (LzZl 

cr... ,  cr. . , 


F2  =  tanh(d>2 ) 


O,  =  max  2 


4k  500/y 
0.09 coy  ’  py2co 


Now  returning  to  the  cross-diffusion  term,  as  the  SST  k-co  model  is  based  upon  both  the 
Standard  k-co  and  the  Standard  k-£  model,  die  k-£  model  has  been  transformed  into 
equations  based  on  k  and  co,  hence  die  cross-diffusion  term,  D„,  This  term  is  defined  as: 


A.  =2(1 

CO  OXj  OXj 
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Unlike  the  Standard  k-co  model,  however,  the  turbulent  viscosity  formulation  is  given 
for  the  high-Reynolds  number  form  of  the  SST  k-co  models  as: 


pk  1 


CO 

max 

1, 

axoo 

(25) 


Where: 

The  constants  applied  in  the  high  Reynolds  number  form  of  the  SST  k-co  turbulence 
model  are  equal  to: 

ok,i  =  1.176,  <rni  =  2.0,  ok, 2  =  1.0,  Oi0l2  =  1.168,  ai  =  0.31,  fk,i—  0.075,  /?,2  =  0.0828,  k—  0.41 


4.  Numerical  Methods 

4.1  Computational  Domain  and  Mesh 

The  computational  domain  extended  1.5C  upstream  of  the  leading  edge  of  the 
hydrofoil,  1.5C  above  and  below  the  pressure  surface  and  3C  downstream  from  the 
trailing  edge.  Velocity  inlet  boundary  conditions  were  applied  upstream,  above  and 
below  the  hydrofoil,  while  downstream  an  outflow  boundary  condition  was  used.  The 
pressure  and  suction  surfaces  of  the  hydrofoil  were  defined  independently  with  no-slip 
wall  boundary  conditions. 

The  selection  of  the  boundary  conditions  and  dimensions  of  the  computational  domain 
were  an  approximation  to  the  free  field.  It  was  concluded  that  the  flow  boundary 
conditions  were  sufficiently  distanced  from  the  hydrofoil's  surfaces  such  they  should 
have  little  impact  on  the  solution  behaviour.  The  effect  of  the  walls  of  the  test  section 
was  considered  by  the  method  discussed  in  [1].  A  free-stream  velocity  correction  was 
calculated  to  be  1.045  for  the  3m/ s  and  1.044  for  6  m/s 

A  finite  volume  method  was  employed  to  obtain  a  solution  of  the  spatially  averaged 
Navier-Stokes  equations.  The  coupling  between  the  pressure  and  velocity  fields  was 
achieved  using  the  SIMPLE  technique.  A  second-order  upwind  scheme  was  used  for 
the  convection  and  the  central-differencing  scheme  for  diffusion  terms. 

An  unstructured  mesh  arrangement  with  quadrilateral  elements  was  adopted  to  map 
the  flow  domain.  Particular  attention  was  directed  to  an  offset  'inner  region' 
encompassing  the  hydrofoil.  Within  this  region  a  considerably  fine  O-type  mesh  was 
applied  to  achieve  sufficient  resolution  of  the  hydrofoil  surface  and  the  boundary  layer 
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region.  Continuing  downstream  from  the  trailing  edge  of  the  inner  region  was  the 
'wake  region'.  A  considerably  fine  H-type  mesh  was  applied  within  this  region  to 
accurately  resolve  the  near  and  far-wake  flow  behaviour.  The  remaining  'outer  region' 
of  the  domain  was  discretized  with  a  much  coarser  H-type  mesh. 

It  should  be  noted  that  on  the  experimental  hydrofoil  used  by  Bourgoyne  et.  al.  [1, 14  & 
15],  the  trailing  edge  was  ground  to  a  0.397mm  radius,  however,  to  allow  improved 
surface  discretization  and  mesh  quality  at  the  trailing  edge  of  the  numerical  model  this 
radius  was  increased  to  1.5mm. 


4.2  Model  Parameters 

The  experimental  hydrofoil  was  CNC  machined  from  a  solid  casting  of  Ni-Al  Bronze. 
The  RMS  surface  roughness  height,  after  polishing,  was  determined  to  be  2.5xl0-7m 
[14],  It  was  assumed  that  the  roughness  distribution  was  'uniform  sand-grain', 
allowing  the  use  of  a  roughness  constant  of  0.5. 

Information  on  the  flow  condition  within  the  LCC's  test  section  during  the 
experimental  testing  was  not  explicitly  detailed  by  Bourgoyne  et.  al.  [1,  14  &  15],  only 
that  it  was  "a  low  turbulence  water  tunnel"  [15].  From  the  literature  [20]  it  was 
determined  that  the  LCC's  Tow  turbulence'  intensity  was  in  the  order  of  0.1%.  This 
value  was  assessed  as  reasonable  and  applied  to  all  applicable  boundary  conditions. 


4.3  Wall  Treatments 

4.3.1  Wall  Functions 

The  application  of  wall  functions  to  modelling  the  near-wall  region  may  significantly 
reduce  both  the  processing  and  storage  requirements  of  a  numerical  model,  while 
producing  an  acceptable  degree  of  accuracy.  The  principal  of  the  wall  function  is  to 
'bridge'  the  viscous  flow  near  a  wall  and  link  the  solution  variables  at  the  wall-adjacent 
cells  to  the  corresponding  quantities  at  the  wall.  With  careful  design  of  near  wall  mesh, 
wall  functions  should  be  applied  where  there  are  several  cells  within  the  boundary 
layer.  Their  use  with  excessively  fine  mesh  near  walls  should  be  avoided  as  their 
validity  expires  within  die  viscous  sublayer  [13], 

With  respect  to  the  non-dimensional  wall  parameter,  i /+,  a  target  value  of  30  is  desired 
of  a  mesh  for  die  application  of  wall  functions.  This  parameter  is  defined  as: 


(26) 

A 


14 


DSTO-TR-1564 


In  Equation  (26),  yp  is  the  distance  from  the  first  computational  node  to  the  wall  and  the 
subscript  w  denotes  wall  properties. 


4.3.1. 1  Standard  Wall  Function 

The  standard  wall  function  applied  is  the  logarithmic  law-of-the-wall  for  mean 
velocity,  given  in  the  literature  [13]  as: 

tr  =-ln(i?/)  (27) 

X* 


Where: 


.  pcT5k°/yP 
y  = - 

P 

kP  is  the  turbulence  kinetic  energy  the  cell  centroidal  node  P  and  yv  is  the  distance  from 
node  P  to  the  wall,  k  and  E  are  constants  equal  to  0.42  and  9.81  respectively. 

For  explicit  detail  on  the  computation  of  turbulence  quantities  in  wall-adjacent  cells 
using  the  standard  wall  function,  the  reader  is  directed  to  the  literature  [13].  Standard 
wall  functions  work  reasonably  well  for  a  wide  range  of  wall  bounded  flows,  however, 
when  the  flow  begins  to  diverge  from  ideal  conditions,  such  as  with  strong  adverse 
pressure  gradients,  the  accuracy  of  their  predictions  is  likely  to  degrade  [13]. 


4.3.1.2  Non-Equilibrium  Wall  Function 

The  main  differences  between  tire  standard  wall  function  and  the  non-equilibrium  wall 
function  are  that  the  law-of-the-wall  for  mean  velocity  is  sensitised  to  pressure  gradient 
effects  and  a  two-layer-based  concept  is  applied  to  calculate  the  cell-averaged 
turbulence  kinetic  energy  production  and  dissipation  in  wall-adjacent  cells  [13]. 

The  law-of-the-wall  for  mean  velocity  sensitised  to  pressure  gradients  is  given  in  the 
literature  [13]  as: 
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yY  is  a  constant  equal  to  11.225.  For  the  cell-averaged  turbulence  kinetic  energy 
production  and  dissipation  in  wall-adjacent  cells,  the  two-layer  concept  assumes  that 
the  wall-adjacent  cells  consist  of  both  a  viscous  sublayer  and  a  fully  turbulent  layer 
[13].  Again  for  explicit  detail  on  the  computation  of  turbulence  quantities  in  these  cells 
using  the  non-equilibrium  wall  function,  the  reader  is  directed  to  the  literature  [13], 


4.3.2  Enhanced  Wall  Treatment 

Enhanced  wall  treatment  is  a  method  of  near-wall  modelling  that  utilises  the 
combination  of  a  two-layer  zonal  model  with  enhanced  wall  functions.  This  wall 
treatment  is  suitable  for  both  considerably  fine  meshes  and  coarser  meshes  (or  'wall 
function  meshes'). 

If  the  mesh  is  fine  enough  to  resolve  the  laminar  sublayer  (within  die  order  of  y+  «  1), 
then  the  wall  treatment  is  identical  to  the  two-layer  zonal  model,  however,  this  mesh 
requirement  can  place  significant  demands  on  computational  processing  and  storage 
infrastructure.  For  full  details  on  the  two-layer  zonal  model,  the  enhanced  wall 
functions  and  their  method  of  blending,  the  reader  is  directed  to  the  literature  [13]. 
Enhanced  wall  treatment  with  near-wall  modelling  is  particularly  required  for  flow 
scenarios  involving  severe  pressure  gradients  leading  to  boundary  layer  separation 


5.  Results  and  Discussion 

5.1  Grid  Independence  Analysis 

A  grid  independence  analysis  was  conducted  using  four  meshes  of  varying  cell 
number.  Each  mesh  was  processed  using  the  Realisable  k-e  turbulence  model  with 
enhanced  wall  treatment,  at  a  free-stream  velocity  of  3m/ s  (Re  =  8.284 xlO6).  Table  1 
details  the  node  and  i p  characteristics  of  each  mesh,  highlighting  the  refinements 
between  the  subsequent  meshes  A,  B,  C  and  D. 

Both  mesh  C  and  mesh  D  achieve  considerably  low  average  i p  values,  sufficiently 
resolving  the  laminar  sublayer  (i.e.  i \p  <  4-5).  Although  having  a  minimum  i p  value 
well  within  the  laminar  sublayer,  the  maximum  and  average  values  of  mesh  B  indicate 
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that  its  resolution  extends  to  the  buffer  layer  only  (i.e.  5  <  \f  <  30).  Also  with  a 
sufficiently  low  minimum  i/+,  mesh  A  shows  significant  coarseness  with  considerably 
larger  maximum  and  average  y+  values,  indicating  resolution  to  the  turbulent  outer 
layer  only  (i.e.  i/+  >  30). 

Table  2  shows  the  predicted  pressure-derived  lift  and  drag  coefficients,  and  the 
respective  percentage  error  from  the  experimental  values  for  each  mesh.  The 
application  of  each  mesh  generally  produces  accurate  predictions  of  the  lift  and  drag 
coefficients  with  comparison  to  the  experimental  data.  It  can  be  seen  that  the  error 
associated  with  the  predicted  lift  coefficient  decreases  with  mesh  refinement,  however, 
at  the  same  time  an  increase  in  error  of  the  predicted  drag  coefficient  is  shown.  With 
the  use  of  mesh  C  and  mesh  D  predicting  almost  identical  coefficients,  this  suggests 
that  the  solution  has  converged  and  grid  independence  has  been  achieved. 

Figure  3  displays  the  pressure  surface  boundary  layer  normalised  mean  velocity 
profiles  at  93%C,  predicted  using  each  mesh.  With  comparison  to  the  experimental 
data,  the  velocity  profile  was  slightly  under-predicted,  but  generally  well  resolved 
using  each  mesh.  The  velocity  profiles  predicted  using  mesh  C  and  mesh  D  are 
identical  and  show  clear  differences  from  those  of  mesh  A  and  mesh  B,  again 
suggesting  grid  independence  with  the  two  finest  meshes. 

Therefore,  with  mesh  C  and  mesh  D  predicting  almost  identical  results,  it  is  concluded 
that  a  suitable  degree  of  grid  independence  is  shown  by  the  numerical  solution.  In 
consideration  of  the  processing  infrastructure  and  time  constraints,  the  very  marginal 
performance  advantages  provided  by  the  finest  mesh  (mesh  D)  do  not  justify  the 
increased  cell  numbers  (~  16%)  and  extended  CPU  processing  times.  Henceforth,  all 
numerical  solutions  discussed  are  obtained  using  mesh  C. 

5.2  Turbulence  Model  Performance 

The  results  of  the  performance  of  each  turbulence  model  are  presented  with  respect  to 
the  predicted  pressure  coefficient  distribution,  lift  and  drag  coefficients,  the  pressure 
surface  boundary  layer,  suction  surface  boundary  layer  separation  and  the  trailing 
edge  and  far-wake  flow.  The  performance  of  each  turbulence  model  is  assessed  at 
chord-based  Reynolds  numbers  of  8.284X106  (Ux=  3m/ s)  and  1.657xl07  (17*.=  6m/ s). 

NOTE:  All  presented  data  is  dimensionless:  co-ordinate  locations  are  referenced  to  the 
chord  length  (i.e.  x/C  and  y/C )  and  all  mean  velocities  are  normalised  with  the  free- 
stream  velocity. 


5.2.1  Coefficients  of  Pressure,  Lift  and  Drag 

Figure  4  illustrates  the  pressure  coefficient  distribution  at  the  surface  of  the  hydrofoil 
as  predicted  by  each  turbulence  model,  at  a  free-stream  velocity  of  3m/  s  (A  negative 
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scale  is  used  to  enable  easier  comparison  with  die  physical  surfaces  of  the  hydrofoil). 
The  distributions  predicted  by  the  Realisable  k-e  and  the  SST  k-co  models  are  very 
encouraging,  however,  the  Standard  k-e  and  Standard  k-co  models  show  significant 
deviation  from  the  experimental  data. 

The  Realisable  k-e  turbulence  model  predicts  high-pressure  (referenced  to  the  free- 
stream  static  pressure)  below  the  hydrofoil  and  low-pressure  above,  whereas  the  SST 
k-co  model  predicts  low-pressure  both  above  and  below.  This  indicates  that  the 
Realisable  k-e  model  predicts  decreased  flow  velocity  below  the  pressure  surface  and 
increased  flow  velocity  above  the  suction  surface,  whereas  the  SST  k-co  predicts 
increased  flow  velocity  both  above  and  below  the  hydrofoil. 

Both  of  these  flow  scenarios  are  valid,  being  dependent  on  the  camber  and  thickness 
characteristics,  however,  as  these  variables  are  fixed,  all  turbulence  models  should 
predict  the  same  flow  behaviour.  Comparison  to  the  experimental  data  demonstrates 
that  the  Realisable  k-e  model  excellently  predicts  the  flow  behaviour  below  the 
pressure  surface.  This  model  also  performs  exceptionally  well  at  the  leading  edge, 
accurately  predicting  the  stagnation  conditions.  However  poorer  correlation  is  shown 
at  the  trailing  edge.  Over  the  larger  suction  surface  the  Realisable  k-e  model  under 
predicts  the  magnitude  of  the  low  static  pressure  region.  The  net  effect  of  this  is  to 
reduce  the  pressure-induced  force  acting  on  the  suction  surface,  hence  the  net  lift  force 
and  the  lift  coefficient. 

The  SST  k-co  turbulence  model  shows  very  similar  qualitative  characteristics  as  to  those 
predicted  by  the  Realisable  k-e  model.  However,  the  predicted  pressure  coefficient 
distribution  is  translated  upwards.  This  translation  results  in  the  SST  k-co  model  more 
closely  predicting  the  pressure  coefficient  distribution  on  the  suction  surface.  However, 
a  poorer  correlation  is  shown  with  the  experimental  data  in  the  leading  and  trailing 
edge  regions,  and  at  the  pressure  surface. 

Similarly,  the  Standard  k-e  and  Standard  k-co  models  predict  almost  identical  qualitative 
characteristics  of  the  pressure  coefficient  distribution  as  to  those  predicted  by  the 
previously  discussed  models.  However,  both  models  significantly  over-predict  the 
pressure  coefficient  distribution,  indicating  vastly  different  behaviour  of  the  flow 
around  the  hydrofoil  with  respect  to  the  experimental  data. 

These  two  models  predict  considerably  low-pressure  regions  both  above  and  below  the 
hydrofoil.  The  results  are  significantly  lower  in  magnitude  than  those  predicted  by 
either  the  Realisable  k-e  or  SST  k-co  turbulence  models.  This  suggests  that  both  models 
predict  significant  flow  accelerations  both  above  and  below  the  hydrofoil.  It  is  expected 
that  this  is  caused  by  the  excessive  level  of  turbulent  viscosity,  which  has  also  been 
observed  by  other  researchers  [4,7], 
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With  such  deviations  from  the  experimental  data,  the  pressure  coefficient  distributions 
predicted  by  the  Standard  k-e  and  Standard  k-oo  turbulence  models  are  viewed  with 
considerable  concern. 

Inspection  of  the  predicted  pressure  coefficient  distribution  at  a  free-stream  velocity  of 
6m/s  (illustrated  in  Figure  5)  reveals  identical  predictions  by  each  turbulence  model  as 
at  3m/ s.  Again  the  Realisable  k-e  turbulence  model  most  accurately  predicts  the 
experimental  data,  while  the  Standard  k-e  and  Standard  k-co  models  show  considerable 
error  in  their  predictions. 

It  is  concluded  that  the  errors  associated  with  the  pressure  coefficient  distributions 
predicted  by  the  Standard  k-e  and  Standard  k-co  turbulence  models  are  not  associated 
with  post-processing  error  for  two  reasons.  Firstly,  the  stagnation  conditions  are 
correctly  predicted  (i.e.  -Cp  «  -1)  by  each  model  and  secondly,  variations  of  similar 
magnitudes  are  present  in  the  reported  pressure  forces  acting  on  the  individual 
surfaces  of  the  hydrofoil,  detailed  in  Table  3  arid  Table  4,  for  free-stream  velocities  of 
3m/ s  and  6m/ s  respectively. 

Although  the  net  lift  forces  predicted  by  each  turbulence  model  are  very  similar,  the 
forces  on  the  individual  surfaces  vary  significantly.  The  Realisable  k-e  model  predicts 
12.36kN  and  49.47kN  acting  on  die  suction  surface  at  free-stream  velocities  of  3m/ s 
and  6m/s  respectively,  however,  the  Standard  k-co  model  predicts  respective  values  of 
51.13kN  and  186.74kN,  values  approximately  four  times  greater.  Again  this 
characteristic  is  in  accordance  with  other  researchers  [4,7]  observations  the  Standard  k-e 
and  Standard  k-co  turbulence  models  over-predict  the  level  of  turbulent  viscosity. 

Therefore,  it  is  assumed  that  diis  is  an  accurate  representation  of  the  performance  of 
each  turbulence  model  in  predicting  the  pressure  coefficient  distribution  of  the 
hydrofoil. 

Table  5  displays  the  pressure-derived  lift  and  drag  coefficients  predicted  by  each 
turbulence  model  at  a  free-stream  velocity  of  3m/ s.  Percentage  error  from  the 
experimental  values  are  indicated.  Each  turbulence  model  generally  acceptably 
predicts  the  lift  coefficient,  although  the  accuracy  trend  of  each  model  shown  in  the 
pressure  coefficient  distribution  was  not  replicated  as  expected.  The  most  accurate 
turbulence  model  was  the  Standard  k-e,  which  predicted  the  hydrofoils  lift  coefficient 
to  within  1.34%  of  the  experimental  value.  As  expected  from  the  analysis  of  the 
pressure  coefficient  distribution  (Figure  4),  the  Realisable  k-e  model  under-predicts  the 
lift  coefficient  by  4.11  % . 

With  consideration  to  the  pressure  coefficient  distribution,  a  better  indication  of  each 
turbulence  model's  performance  is  gained  from  the  predicted  drag  coefficients.  It  is 
evident  that  the  errors  associated  with  the  Standard  k-e  and  Standard  k-co  turbulence 
models  in  predicting  the  pressure  coefficient  distribution  (although  not  greatly 
affecting  the  predicted  coefficients  of  lift)  have  considerable  impact  on  the  predicted 
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drag  coefficients.  The  accuracy  trend  of  the  turbulence  models  is  identical  to  that 
shown  by  the  pressure  coefficient  distribution;  the  Realisable  k-e  model  predicts  the 
drag  coefficient  most  accurately  (to  within  11.11%  of  the  experimental  value)  and  the 
Standard  k- a  model  the  least  accurately  (444.44%). 

Table  6  displays  the  pressure-derived  lift  and  drag  coefficients  predicted  by  each 
turbulence  model  at  a  free-stream  velocity  of  6m/s.  At  this  flow  velocity,  each 
turbulence  model  reasonably  predicts  the  lift  coefficient  and  again  a  misleading 
performance  trend  is  shown. 

Comparison  to  the  numerical  data  at  a  free-stream  velocity  of  3m/s  (Table  5)  shows 
that  an  increase  in  Reynolds  number  results  in  improvements  in  the  performance  of 
each  turbulence  model.  The  magnitude  of  error  associated  with  the  lift  coefficient 
predicted  by  the  Standard  k-e,  the  Realisable  k-e,  the  Standard  k-co  and  the  SST  k-co 
models  decreased  by  0.31%,  2.63%,  2.8%  and  2.28%  respectively.  For  the  predicted  drag 
coefficient,  however,  the  error  associated  with  both  the  Standard  k-e  and  SST  k-co 
models  increased  by  13.02%  and  5.69%  respectively,  whereas  that  of  the  Realisable  k-e 
and  Standard  k-co  models  showed  slightly  improved  predictions  by  0.26%  (11.11%  to 
10.85%)  and  40.58%  (444.44%  to  403.87%)  respectively. 

In  summary,  the  presented  data  of  the  predicted  pressure  coefficient  distribution  and 
the  lift  and  drag  coefficients  at  free-stream  velocities  of  3m/s  (Re  =  8.284x106)  and 
6m/s  (Re  =  1.657xl07),  the  Realisable  k-e  turbulence  model  clearly  outperforms  the 
remaining  turbulence  models.  This  model  also  consistently  shows  equal,  if  not  slightly 
improved,  predictions  at  the  higher  Reynolds  number. 


5.2.2  Pressure  Surface  Boundary  Layer 

Figure  6  illustrates  die  numerical  and  experimental  pressure  surface  boundary  layer 
velocity  profiles  at  93%C  at  a  free-stream  velocity  of  3m/s.  The  results  have  been 
inverted  to  enable  comparison.  As  the  pressure  surface  is  largely  flat  beyond  28%C 
with  only  moderate  surface  gradients  prior  to  this  point,  it  can  be  approximated  as  a 
flat  plate.  Therefore,  also  shown  in  Figure  6  is  the  theoretically  predicted  turbulent 
boundary  layer  velocity  profile,  calculated  using  classical  turbulent  boundary  layer 
equations  [21]. 

Each  turbulence  model  predicts  a  differing  velocity  profile  widi  varying  degrees  of 
correlation  to  the  experimental  data.  Good  correlation  is  shown  by  the  theoretical 
velocity  profile  witii  accurate  predictions  of  botii  the  velocity  gradient  and  the 
boundary  layer  thickness  (defined  as  the  endpoint  of  the  profile). 

The  Standard  k-e,  Standard  k-co  and  SST  k-co  turbulence  models  all  produce  smoothly 
changing  profiles  that  slightly  over-predict  the  velocities  closer  to  the  surface  (i.e.  y/C 
<  0.009)  and  under-predict  further  away  from  the  surface. 
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The  turbulence  model  that  most  accurately  predicts  the  boundary  layer  velocity  profile 
is  the  Realisable  k-s  model.  Close  to  the  surface  (i.e.  y/C  <  0.005)  this  model  predicts 
the  boundary  layer  characteristics  exceptionally  well  and  continues  to  maintain  a  high 
degree  of  correlation  further  from  the  surface.  Unlike  the  other  turbulence  models,  the 
Realisable  k-s  model  produces  a  definable  gradient  change  where  the  velocity  profile 
becomes  quite  blunt.  The  approximate  location  where  the  velocity  profile  exhibits  this 
characteristic  visually  correlates  well  to  the  experimental  and  theoretical  boundary 
layer  thicknesses  (i.e.  u  «  99%UX). 

Table  7  details  the  boundary  layer  parameters  derived  from  the  numerical  data  of  each 
turbulence  model,  and  presents  a  comparison  with  the  experimental  and  theoretical 
parameters  at  a  free-stream  velocity  of  3m/ s. 

The  superior  performance  of  the  Realisable  k-e  model  is  highlighted  in  Table  7.  The 
boundary  layer  thickness  and  displacement  thickness  are  predicted  within  6.18%  and 
9.60%  respectively.  The  SST  k-co  model  predicts  these  parameters  to  within  27.68%  and 
16.16%  respectively,  whilst  the  Standard  k-e  and  Standard  k-co  turbulence  models 
significantly  over-predicted  the  boundary  layer  thickness  (164.25%  and  291.02% 
respectively)  and  over-predict  the  displacement  thickness  (16.86%  and  15.69% 
respectively). 

Figure  7  illustrates  the  same  pressure  surface  boundary  layer  profiles  at  93  %C,  at  a 
free-stream  velocity  of  6m/ s.  Again  the  Realisable  k-e  turbulence  model  excellently 
predicts  the  boundary  layer  velocity  profile.  Slight  over-predictions  are  noted  closer  to 
the  surface  (i.e.  y/C  <  0.009).  A  poorer  correlation  to  the  experimental  data  is  shown  by 
the  theoretically  calculated  velocity  profile,  as  the  Reynolds  number  of  the  flow  is 
slightly  greater  than  the  valid  Reynolds  number  range  of  the  equation  [21].  The 
Standard  k-s,  Standard  k-co  and  SST  k-co  turbulence  models  again  all  produce  smoothly 
changing  velocity  profiles  that  over-predict  the  experimental  velocity  profile  closer  to 
the  surface. 

Similarly  to  Table  7,  Table  8  details  the  boundary  layer  parameters  at  a  free-stream 
velocity  of  6m/s.  Given  that  the  uncertainty  of  the  experimental  data,  measured  a 
boundary  thickness  for  6m/s  greater  than  that  for  3m/s,  the  comparison  with 
experimental  data  seems  unnecessary.  However,  in  Table  8,  the  experimental  data  was 
listed  as  a  reference.  All  models  correctly  predict  decreasing  boundary  layer  thickness 
with  increasing  Reynolds  number,  except  the  SST  k-co  turbulence  model,  which 
predicts  a  greater  boundary  layer  thickness  at  6m/ s  (38.06mm)  compared  to  that  at 
3m/ s  (37.82mm). 

Both  the  Standard  k-e  and  Standard  k-co  turbulence  models  again  performed 
considerably  poorly,  significantly  over-predicting  both  the  boundary  layer  thicknesses 
and  the  displacement  thicknesses.  Like  the  SST  k-co  model  these  models  incorrectly 
predict  the  trend  of  the  boundary  layer  thickness  with  increasing  Reynolds  number. 
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In  summary  of  the  predictions  of  the  pressure  surface  boundary  layer  at  93  %C  at  free- 
stream  velocities  of  both  3m/ s  and  6m/s,  the  Realizable  k-e  turbulence  model  clearly 
outperforms  the  Standard  k-e,  Standard  k-co  and  SST  k-co  models.  This  model  predicts 
the  velocity  profile,  boundary  layer  thickness  and  displacement  thickness  to  a  high 
degree,  and  was  the  sole  turbulence  model  to  predict  boundary  layer  thinning  with 
increasing  Reynolds  number. 


5.2.3  Suction  Surface  Boundary  Layer  Separation 

To  assess  if  each  turbulence  model  correctly  predicts  suction  surface  boundary  layer 
separation,  analysis  of  the  coefficient  of  skin  friction  on  the  surfaces  of  the  hydrofoil 
was  completed.  For  a  two-dimensional  steady-state  flow,  the  full  separation  point  (S) 
of  the  boundary  layer  can  be  approximately  defined  where  the  wall  shear  stress,  hence 
the  skin-friction  coefficient,  becomes  zero  [21], 

Figure  8  illustrates  the  skin-friction  coefficient  distribution  on  the  suction  surface  at  the 
trailing  edge  for  a  free-stream  velocity  of  3m/s.  Experimentally,  boundary  layer 
separation  was  determined  to  occur  between  95.75%C  and  97.50%C  on  the  suction 
surface  of  the  hydrofoil  [1], 

Both  the  Realisable  k-e  and  SST  k-co  turbulence  models,  and  the  Standard  k-e  and 
Standard  k-co  turbulence  models  predict  different  skin-friction  coefficient  magnitudes 
in  the  vicinity  of  the  trailing  edge.  However,  distinctive  features  predicted  by  each 
turbulence  model  occur  at  very  similar  chord-wise  locations. 

At  approximately  96.5%C,  all  models  predict  a  fluctuating  increase  in  the  skin-friction 
coefficient,  peaking  at  a  chord-wise  location  of  96.82%C.  It  is  considered  that  these 
fluctuations  are  associated  with  the  combined  effects  of  the  considerably  fine  mesh  and 
the  surface  discontinuities  introduced  when  the  trailing  edge  was  modified.  Directly 
following  the  skin-friction  coefficient  peak  is  a  steep  decrease,  reaching  zero  at  a  chord- 
wise  location  of  between  97.86%C  and  98.52%C  (depending  on  the  turbulence  model). 
Beyond  this  chord-wise  location,  a  small  secondary  increase  of  the  skin-friction 
coefficient  occurs,  followed  by  a  gradual  decrease  to  zero  at  the  trailing  edge  (In  Figure 
8,  the  trailing  edge  terminates  just  prior  to  100%C  due  to  the  fillet  radius). 

It  is  considered  that  the  fluctuating  increase  then  sharp  decrease  in  the  skin  friction 
coefficient  is  associated  with  the  commencement  of  suction  surface  boundary  layer 
separation  at  die  commencement  of  die  trailing  edge  modifications  at  96.75%C.  The 
steep  adverse  pressure  gradient  widiin  this  region  begins  to  significantly  retard  the 
motion  of  fluid  particles  adjacent  to  the  surface  reducing  their  velocity,  hence  the  wall 
shear  stress  and  the  coefficient  of  skin  friction. 
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The  first  chord-wise  location  where  the  skin  friction  coefficient  steeply  decreases  to 
zero  is  considered  to  be  the  location  at  which  the  boundary  layer  flow  adjacent  to  the 
surface  stagnates.  This  indicates  that  the  fluid  particles  within  the  boundary  layer  flow 
directly  adjacent  to  the  surface  are  halted  through  the  combined  effects  of  the  adverse 
pressure  gradient  and  backflow  from  the  recirculation  bubble.  The  second  chord-wise 
location  where  the  skin-friction  coefficient  tends  to  zero  again  indicates  flow  stagnation 
and  the  partial  reattachment  of  the  boundary  layer  flow.  The  chord-wise  distance 
between  these  two  points  of  boundary  layer  flow  stagnation  defines  the  length  of  the 
recirculation  bubble. 

As  all  turbulence  models  predict  the  commencement  of  suction  surface  boundary  layer 
separation  at  the  start  of  the  modified  trailing  edge  (i.e.  «96.75%C),  the  only  variation 
between  each  models  performance  is  the  prediction  of  the  chord-wise  location  of  full 
boundary  layer  separation  (i.e.  the  first  location  of  zero  skin-friction  and  boundary 
layer  flow  stagnation),  hence  the  chord-wise  size  of  the  recirculation  bubble. 

The  Standard  k-e,  Realisable  k-e,  Standard  k-co  and  SST  k-co  turbulence  models  predict 
full  boundary  layer  separation  at  chord-wise  locations  of  98.59%C,  98.37%C,  98.31  %C 
and  98.15%C  respectively.  Therefore,  the  Standard  k-e  model  predicts  the  smallest 
recirculation  bubble,  whereas  the  SST  k-co  predicts  the  largest.  Although  these  locations 
are  slightly  over-predicted  with  comparison  to  the  separation  location  range  predicted 
experimentally  (i.e.  95.75%C  <  S  >  97.50%C),  the  point  at  which  boundary  layer 
separation  commences  correlates  well. 

Figure  9  illustrates  the  skin-friction  coefficient  distribution  on  the  suction  surface  at  tire 
trailing  edge  for  a  free-stream  velocity  of  6m/  s.  Similar  to  the  coefficient  of  skin  friction 
characteristics  at  3m/ s,  the  Realisable  k-e  and  SST  k-co  turbulence  models,  and  the 
Standard  k-e  and  Standard  k-co  turbulence  models  again  predict  different  magnitudes  in 
the  vicinity  of  the  trailing  edge.  Again  the  rapid  decrease  in  the  skin-friction  coefficient, 
indicating  the  commencement  of  boundary  layer  separation,  is  predicted  at 
approximately  96.55%C  by  each  turbulence  model.  This  suggests  that  the  separation 
commencement  location  is  independent  of  Reynolds  number  and  that  the  trailing  edge 
geometry  is  the  dominant  variable  in  its  location. 

Similarly,  two  locations  where  the  skin  friction  coefficient  equalled  zero  were  predicted 
on  the  trailing  edge  curve.  The  predicted  locations  of  full  boundary  layer  separation  by 
the  Standard  k-e,  Realisable  k-e,  Standard  k-co  and  SST  k-co  turbulence  models,  are 
98.75%C,  98.48%C,  98.48%C  and  98.31%C  respectively.  When  compared  with  the 
respective  chord- wise  locations  at  3m/ s,  each  location  of  full  boundary  layer 
separation  shifts  closer  to  the  trailing  edge  by  0.16%C,  0.11  %C,  0.17%C  and  0.16%C 
respectively.  Again  the  Standard  k-e  predicted  the  smallest  recirculation  bubble  and  the 
SST  k-co  model  predicted  the  largest. 

This  suggests  that  with  increasing  Reynolds  number  and  a  fixed  commencement  point 
of  boundary  layer  separation,  the  location  at  which  the  boundary  layer  flow  fully 
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separates  moves  rearwards,  confirming  the  experimental  observations  of  both 
Bourgoyne  et.  al.  [1]  and  Baubeau  &  Latorre  [17]. 

In  summary,  each  turbulence  model  accurately  predicts  the  commencement  location  of 
boundary  layer  separation  and  that  the  location  of  full  separation  moves  rearwards 
with  increasing  Reynolds  number.  Although  it  is  not  possible  to  determine  the 
turbulence  model  that  provides  superior  performance  with  respect  to  the  experimental 
data,  the  solution  convergence  shown  by  both  the  Realisable  k-e  and  SST  k-co  models  at 
both  3m/s  and  6m/s  indicates  a  more  confident  prediction  of  the  skin  friction 
coefficient  in  the  trailing  edge  region. 


5.2.4  Trailing  Edge  and  Wake  Flow 

To  assess  the  performance  of  each  turbulence  model  in  predicting  the  flow 
characteristics  in  the  near  and  far  wake,  the  normalised  stream-wise  mean  velocity 
profiles  at  three  chord-wise  planes  in  the  near  wake  and  one  in  the  far  wake  are 
discussed  for  free-stream  velocities  of  3m/s  and  6m/s. 

Figure  10  illustrates  the  normalised  stream-wise  mean  velocity  profile  within  the 
trailing  edge  region  at  chord-wise  location  of  97.8%C.  This  plane  is  located  between  the 
start  of  the  trailing  edge  modifications  (96.75%C)  and  the  tip  of  the  trailing  edge,  prior 
to  the  point  of  full  boundary  layer  separation. 

The  Standard  k-co,  Standard  k-e,  and  the  SST  k-co  models  generally  perform  poorly 
within  this  region,  under-predicting  the  velocities  above  and  below  the  trailing  edge 
surfaces,  possibly  suggesting  that  the  models  are  too  diffusive.  The  Realisable  k-e 
model,  however,  shows  encouraging  performance  in  resolving  the  important  flow 
features.  This  model  predicts  the  steep  velocity  gradient  changes  through  the 
separating  boundary  layer  flow.  Below  the  pressure  surface,  improved  performance  is 
shown  by  the  Realisable  k-e  model  in  predicting  the  pressure  surface  boundary  layer 
flow  near  the  trailing  edge  tip.  This  model  is  the  only  model  to  accurately  predict  a 
slight  'surge'  in  the  velocity  profile,  located  near  y/C»  0.018. 

Figure  11  illustrates  the  normalised  stream-wise  mean  velocity  profile  within  the 
trailing  edge  region  at  die  tip  of  the  hydrofoil  (i.e.  100%C),  a  region  of  considerable 
turbulent  flow  activity.  Although  accurately  predicting  the  magnitude  of  the  flow 
reversal  associated  with  the  recirculation  bubble,  the  Standard  k-e,  Standard  k-co  and 
SST  k-co  turbulence  models  show  poor  correlation  to  the  experimental  data.  Again 
neither  of  these  models  predict  the  steep  velocity  gradient  changes  through  the 
recirculation  bubble  and  below  the  pressure  surface.  The  Standard  k-e  and  Standard  k- 
o)  models  show  considerable  error  above  the  suction  surface. 

The  Realisable  k-e  model  accurately  predicts  the  magnitude  of  the  flow  reversal  and 
shows  excellent  correlation  to  the  experimental  data  below  the  trailing  edge,  again 
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accurately  predicting  the  velocity  profile  'surge'.  Above  the  trailing  edge,  the 
correlation  to  the  experimental  data  degrades  marginally  with  an  over  prediction  in  the 
transverse  thickness  of  the  recirculation  bubble. 

Figure  12  illustrates  the  normalised  stream-wise  mean  velocity  profile  just  downstream 
of  the  hydrofoil's  trailing  edge  at  102.8 %C,  within  the  near  wake  region.  In  general,  a 
poorer  velocity  profile  prediction  is  shown  by  each  turbulence  model  compared  with 
previous  chord-wise  locations.  All  models  over-predict  the  minimum  velocity  at  the 
core  of  the  wake  and  under-predict  the  transverse  thickness  of  the  wake  region, 
showing  poor  correlation  to  the  experimental  data  above  the  trailing  edge.  The 
Realisable  k-s  turbulence  model  is  the  only  model  to  successfully  predict  the  steep 
velocity  gradient  changes  through  the  wake  core.  However,  beyond  the  wake  core  it 
slightly  over  predicts  the  velocity  profile. 

Figure  13  illustrates  the  normalised  stream-wise  mean  velocity  at  150%C  within  the  far 
wake,  a  region  of  continuing  turbulent  kinetic  energy  dissipation.  When  compared 
with  the  turbulence  model  performance  observed  in  the  near  wake  (Figure  12)  and  the 
over-prediction  of  the  minimum  velocity,  all  models  in  the  far  wake  under-predict  the 
minimum  velocity  within  the  wake  core.  Again  the  velocity  profile  predicted  by  the 
Standard  k-s  and  Standard  k-co  models  are  considerably  poor.  The  Realisable  k-s  model 
predicts  almost  instantaneous  velocity  gradient  changes  through  the  wake  core,  while 
the  SST  k-  co  model  predicts  steep,  yet  more  gradual  velocity  gradient  changes, 
suggesting  that  both  models  posses  very  low  damping.  Both  models  generally  show 
good  correlation  to  the  experimental  data,  particularly  within  the  wake  core.  The 
'downwash'  effect  of  the  flow  over  the  hydrofoil  is  generally  well  predicted  by  each 
turbulence  model,  with  the  wake  core  being  translated  downwards  with  respect  to  its 
previous  position  in  the  near  wake  (Figure  12). 

For  analysis  of  the  trailing  edge  and  wake  flow  at  a  higher  Reynolds  number.  Figure 
14,  Figure  15,  Figure  16  and  Figure  17  illustrate  the  normalised  stream-wise  mean 
velocity  profiles  at  97.8%C,  100%C,  102.8%C  and  150%C  respectively,  at  a  free-stream 
velocity  of  6m/s,  correlating  to  Re  =  1.657xl07. 

When  compared  with  the  performance  predicted  at  the  low  velocity  (Figure  10),  a 
distinct  improvement  is  achieved  by  the  Realisable  k-s  turbulence  model  (Figure  14).  It 
can  be  seen  in  Figure  14  that  above  the  suction  surface  the  Realisable  k-s  model  shows 
considerable  improvement  in  predicting  tire  separating  boundary  layer  flow,  while 
below  the  pressure  surface  the  pressure  surface  boundary  layer  is  more  accurately 
resolved.  The  results  presented  in  Figure  10  indicates  that  the  remaining  turbulence 
models  have  considerable  difficulty  in  predicting  the  steep  velocity  gradient  changes 
through  the  boundary  layer  flow  both  above  and  below  the  trailing  edge  surfaces. 

Similar  performance  improvements  are  shown  by  the  Realisable  k-s  model  in  Figure  15, 
at  a  chord-wise  location  of  100  %C.  The  accuracy  of  the  predicted  flow  reversal  in  the 
recirculation  bubble  is  lower  than  at  3m/ s,  however  a  higher  degree  of  correlation  to 
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the  experimental  velocity  profile  is  shown  through  the  recirculation  bubble  and  the 
pressure  surface  boundary  layer.  Although  the  SST  k-co  model  more  accurately  predicts 
the  degree  of  flow  reversal,  it  again  has  difficulty  resolving  the  steep  velocity  gradient 
changes  through  both  the  recirculation  bubble  and  the  pressure  surface  boundary 
layer. 

Within  the  near  wake  at  a  chord-wise  location  of  102.8 %C  (Figure  16),  again  a  generally 
a  poorer  velocity  profile  prediction  is  shown  by  each  turbulence  model  compared  to 
the  previous  chord-wise  planes  (Figure  14  and  Figure  15).  All  models  again  over¬ 
predict  the  minimum  velocity  at  the  core  of  the  wake  and  under-predict  the  transverse 
thickness  of  the  wake  region,  with  the  poorest  correlation  to  the  experimental  data 
occurring  above  the  trailing  edge.  A  considerable  improvement  is  shown  by  the 
Realisable  k-e  model  in  predicting  the  velocity  gradient  changes  through  the  wake  core 
and  the  velocities  beyond  the  wake  core. 

Figure  17  illustrates  the  normalised  stream-wise  mean  velocity  profile  in  the  far  wake 
region  at  150%C.  Again  as  present  within  the  near  and  far  wake  velocity  profiles  at 
3m/ s  (Figures  10, 11, 12,  and  13),  it  appears  that  the  performance  of  each  turbulence 
model  degrades  within  the  wake  region  with  increasing  distance  downstream  from  the 
hydrofoil.  Again  the  Standard  k-e,  Standard  k-co  and  the  SST  k-co  models  show  poor 
correlation  to  the  experimental  velocity  profile.  The  Realisable  k-e  model  acceptably 
predicts  the  experimental  velocity  profile  within  the  wake  core,  with  consideration  to 
the  applied  scale.  Beyond  the  wake  core,  marginal  over-predictions  in  the  velocity 
profile  are  shown. 

In  summary,  when  comparing  the  performance  of  each  turbulence  model  in  predicting 
the  flow  characteristics  in  the  near  and  far  wake,  it  is  apparent  that  the  Realisable  k-e 
turbulence  model  predicts  superior  velocity  profiles.  The  distinct  advantage  of  this 
model  compared  with  the  Standard  k-e,  Standard  k-co  and  SST  k-co  models,  is  that  the 
steep  velocity  gradient  changes  through  the  separating  boundary  layer  flow,  the 
recirculation  bubble,  and  in  the  near  and  far  wakes,  are  accurately  resolved.  Although 
the  Realisable  k-e  model  produces  accurate  predictions  at  3m/s  (Re  =  8.284x106),  it 
shows  considerably  improved  performance  at  6m/s  (Re  =  1.657X107). 

5.3  Analysis  of  Wall  Treatments 

The  results  of  the  performance  of  the  Realisable  k-e  turbulence  model  using  the  three 
wall  treatments  (the  standard  and  non-equilibrium  variants,  and  the  enhanced  wall 
treatment  employing  near-wall  modelling  techniques),  are  presented  at  a  free-stream 
velocity  of  3m/s  (Re  =  8.284x106).  The  results  presented  include  the  predicted  pressure 
coefficient  distribution,  the  lift  and  drag  coefficients,  and  the  pressure  surface 
boundary  layer. 
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Figure  18  illustrates  the  coefficient  of  pressure  distribution  at  tire  surface  of  the 
hydrofoil,  predicted  using  the  diree  different  wall  treatments.  It  is  clear  that  the  use  of 
the  wall  treatments  impacts  upon  the  solution  behaviour.  As  expected  in  flow  scenarios 
involving  strong  adverse  pressure  gradients  and  boundary  layer  separation,  the 
enhanced  wall  treatment  produces  the  most  accurate  distribution.  The  use  of  the 
standard  wall  function  produces  a  slightly  less  accurate  prediction,  yet  performs 
generally  well  at  the  leading  and  trailing  edges.  However,  the  use  of  the  non¬ 
equilibrium  wall  functions  produces  questionable  results,  particularly  at  the  leading 
edge  where  the  coefficient  of  pressure  distributions  on  the  suction  and  pressure 
surfaces  cross  each  other. 

Similarly  to  the  coefficient  of  pressure  distribution,  the  wall  treatment  has  a  marked 
effect  on  the  predicted  lift  and  drag  coefficients.  Table  9  displays  the  pressure-derived 
coefficients  of  lift  and  drag  predicted  by  each  wall  treatment  at  a  free-stream  velocity  of 
3m/ s. 

Using  the  standard  wall  function,  the  coefficient  of  lift  is  largely  under-predicted  by 
17.95%.  However,  this  wall  treatment  produces  a  surprisingly  accurate  prediction  of 
the  coefficient  of  drag,  only  under-predicting  the  experimental  value  by  3.70%.  Much 
poorer  predictions  are  shown  using  the  non-equilibrium  wall  functions,  under¬ 
predicting  the  coefficient  of  lift  by  28.75%  and  over  predicting  the  coefficient  of  drag  by 
33.33%.  The  use  of  enhanced  wall  treatments  produces  generally  accurate  predictions 
of  both  the  lift  and  drag  coefficients  to  within  4.11%  and  11.11%  respectively. 

Figure  19  illustrates  the  numerical,  experimental  and  theoretical  pressure  surface 
boundary  layer  velocity  profiles  at  93  %C,  at  a  free-stream  velocity  of  3m/ s.  As 
expected  in  the  near  wall  region,  the  use  of  each  wall  treatment  predicts  a  different 
boundary  layer  velocity  profile. 

Again  the  use  of  enhanced  wall  treatment  produces  the  most  accurate  velocity  profile 
when  compared  with  the  experimental  data.  The  use  of  the  standard  wall  function 
produces  a  boundary  layer  velocity  profile  with  similar  features  to  those  produced  by 
the  enhanced  wall  treatment,  however  a  larger  boundary  layer  thickness  is  predicted. 
Again,  similar  features  are  predicted  using  the  non-equilibrium  wall  function  and, 
again  a  much  greater  boundary  layer  thickness  is  predicted.  The  use  of  the  enhanced 
wall  treatment  predicts  the  boundary  layer  thickness  to  within  6.18%,  whereas  the 
standard  wall  function  and  non-equilibrium  wall  function  produce  predictions  to 
within  15.26%  and  37.54%  respectively. 

These  results  obtained  using  the  three  different  wall  treatments  reaffirm  the  guidelines 
for  the  application  of  wall  treatments  to  certain  meshes  and  flow  scenarios.  With  an 
average  y+  value  of  2.31,  the  mesh  used  for  this  analysis  is  particularly  fine  (less  than 
the  desired  value  of  30)  and  is  not  classed  as  a  'wall  function  mesh'.  Therefore,  the 
validity  of  wall  functions  applied  to  this  mesh  is  significantly  degraded,  producing 
poor  numerical  results.  It  is  apparent  that  the  standard  wall  function  copes  best  in  such 
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an  adverse  application,  whereas  the  non-equilibrium  wall  function  performs  rather 
poorly. 


6.  Concluding  Remarks 

Evaluation  of  the  Standard  k-s,  Realisable  k-s,  Standard  k-co  and  the  SST  k-co  RANS 
turbulence  models  against  high  quality  experimental  data,  has  revealed  that  the 
Realisable  k-s  model  shows  vastly  superior  performance  when  applied  in  the  steady- 
state  analysis  of  hydrofoil  turbulent  boundary  layer  separation  at  high  Reynolds 
numbers  of  8.284x106  (3m/s)  and  1.657x106  (6m/s). 

The  Realisable  k-£  turbulence  model  accurately  predicts  the  pressure  coefficient 
distribution  at  the  surface  of  the  hydrofoil,  leading  to  good  predictions  in  the  overall 
pressure-derived  lift  and  drag  coefficients.  Analysis  of  the  turbulent  boundary  layer  on 
the  pressure  surface  of  the  hydrofoil  at  93%C  shows  that  the  Realisable  k-s  model 
resolves  the  velocity  profile,  the  boundary  layer  thickness  and  displacement  thickness 
to  a  high  degree.  This  model  is  also  the  only  model  to  correctly  predict  boundary  layer 
thinning  with  an  increase  in  Reynolds  number. 

Suction  surface  boundary  layer  separation  from  the  trailing  edge  is  also  well  predicted 
by  the  Realisable  k-s  turbulence  model.  The  commencement  location  of  boundary  layer 
separation  is  accurately  predicted,  and  this  model  correctly  predicts  that  the  full 
separation  point  of  the  turbulent  boundary  layer  moves  rearward  with  increasing 
Reynolds  number. 

Even  in  the  challenging  separated  trailing  edge  and  near  wake  regions,  the  Realisable 
k-s  model  accurately  predicts  the  normalised  stream-wise  mean  velocity  profiles  at  four 
chord-wise  planar  locations.  The  distinct  performance  advantage  offered  by  the 
Realisable  k-s  model  at  these  locations  that  it  accurately  predicts  the  steep  velocity 
gradient  changes  through  the  separating  boundary  layer  flow,  the  recirculation  bubble 
and  the  wake  core  in  the  near  region. 

The  comparison  of  the  respective  results  at  Reynolds  numbers  of  8.284x106  (3m/s)  and 
1.657x106  (6m/s)  indicates  that  the  performance  of  the  Realisable  k-s  model  shows 
improvement  with  an  increase  in  Reynolds  number.  Predictions  of  tire  hydrofoils 
overall  coefficients  of  lift  and  drag  show  improvement,  the  pressure  surface  boundary 
layer  parameters  are  more  accurately  predicted,  and  considerable  improvement  is 
shown  in  the  prediction  of  the  normalised  stream-wise  mean  velocity  profiles  in  the 
near  and  far  wakes. 
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Furthermore  analysis  of  three  different  wall  treatments  (standard  and  non-equilibrium 
wall  functions,  and  enhanced  wall  treatment  with  near-wall  modelling),  indicates  that 
for  flows  involving  severe  adverse  pressure  gradients  and  boundary  layer  separation. 
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the  application  of  enhanced  wall  treatment  may  be  preferable  for  producing  accurate 
results.  Confirming  the  guidance  provided  in  the  literature  [13],  the  application  of  wall 
functions  to  a  considerably  fine  mesh  significantly  degrades  their  validity,  producing 
poor  numerical  results. 

Therefore,  the  application  of  the  Realisable  k-s  turbulence  model  with  enhanced  wall 
treatment  to  high-Reynolds  number  hydrodynamic  flow  scenarios,  involving  adverse 
pressure  gradients  and  boundary  layer  separation,  may  provide  researchers  with  more 
encouraging  results  from  numerical  methods. 
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NACA-16  Series  Airfoil 
fm  =  3.2%C 


Figure  1.  The  two-dimensional  hydrofoil  geometry  as  used  by  Bonrgoyne  et.  al.  [1,  14  and  15] 
(Adapted  from  1) 
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Figure  2.  The  dimensions  and  boundary  conditions  of  the  computational  domain,  including  the 
dimensions  of  the  water  tunnel's  test  section  (in  which  the  hydrofoil  was  centrally 
positioned) 
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Figure  3.  Grid  Independence  Analysis:  Pressure  surface  boundary  layer  normalised  stream 
mean  velocity  profile  at  93% C  (U^-  3m/s) 
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Figure  4.  Turbulence  Model  Performance: 
the  hydrofoil  (IL  =  3m/ s) 
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coefficient  (Cp)  distribution  at  the  surface  of 
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Figure  5.  Turbulence  Model  Performance:  Pressure  coefficient  (CP)  distribution  at  the  surface  of 
the  hydrofoil  (Ur  =  6m/s) 


Figure  8.  Turbulence  Model  Performance:  Skin-friction  coefficient  distribution  on  the  suction 
surface  at  the  trailing  edge  (Ux  =  3m/s) 
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Figure  9.  Turbulence  Model  Performance:  Skin-friction  coefficient  distribution  on  the  suction 
surface  at  the  trailing  edge  (U,  =  6m/ s) 


Figure  10.  Turbulence  Model  Performance:  Normalised  stream-wise  mean  velocity  profile  at 
97.8% C  (Ux  =  3m/s) 


Ux/Uref 


Figure  17.  Turbulence  Model  Performance:  Normalised  stream-wise  mean  velocity  profile  at 
150%  C  (Ux  =  6m/s) 
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Figure  18.  Wall  Treatment  Analysis:  Pressure  coefficient  (Cp)  distribution  at  the  surface  of  the 
hydrofoil  (Ux  =  3m/s) 
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Table  1.  Grid  Independence  Analysis  -  Mesh  node  and  y+  characteristics 


1 

Nodes 

Surface 

Nodes 

Maximum 

Minimum 

r 

Average 

Mesh  A 

42,282 

570 

220.11 

3.37 

104.35 

MeshB 

98,704 

1060 

31.05 

0.44 

14.30 

MeshC 

208,416 

1590 

4.06 

0.09 

2.31 

MeshD 

244,026 

1590 

1.22 

0.01 

0.68 

Table  2.  Grid  Independence  Analysis  -  Pressure-derived  Lift  and  Drag  coefficients 


Lift  coefficient 

Drag  coefficient 

Experimental 

0.5520 

0.0027 

Predicted 

Error 

Predicted 

Error 

Mesh  A 

0.4945 

10.42% 

0.0029 

7.41% 

MeshB 

0.5230 

5.25% 

0.0025 

7.41% 

MeshC 

0.5302 

3.95% 

0.0024 

11.11% 

MeshD 

0.5305 

3.89% 

0.0024 

11.11% 

Table  3.  Turbulence  Model  Performance:  Hydrofoil  surface  forces  (Ux  =  3m/ s) 


p  .  m/A'-  '''' 

Surface 

Turbulence  Model 

Standard  k- 

£ 

Realizable 

k-£ 

Standard  k~ 

(0 

SST  k-co 

Pressure 

-9.75  kN 

3.09  kN 

-36.23  kN 

-0.24  kN 

Suction 

25.63  kN 

12.36  kN 

51.13  kN 

15.30  kN 

Net  lift  force 

15.88  kN 

15.45  kN 

14.90  kN 

15.54  kN 

Table  4.  Turbulence  Model  Performance:  Hydrofoil  surface  forces  (Ux  =  6m/s ) 


Surface 

Turbulence  Model 

i  -  '  # 

Standard  k- 

£ 

Realizable 

k-£ 

Standard  k- 

(O 

SST  k-w 

Pressure 

-38.66  kN 

13.54  kN 

-125.77  kN 

-0.38  kN 

Suction 

103.28  kN 

49.47  kN 

186.74  kN 

61.55  kN 

Net  lift  force 

64.62  kN 

63.01  kN 

60.97  kN 

61.93  kN 
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Table  5.  Turbulence  Model  Performance:  Pressure-derived  Lift  and  Drag  coefficients 
(Ux  =  3m/s) 


Lift  coefficient 

Drag  coefficient 

Experimental 

0.5520 

0.0027 

Predicted 

Error 

Predicted 

Error 

Standard  k-e 

0.5447 

1.34% 

0.0063 

133.33% 

Realisable  k-e 

0.5302 

4.11% 

0.0024 

11.11% 

Standard  k-o 

0.5113 

7.96% 

0.0147 

444.44% 

SST  k-co 

0.5167 

6.83% 

0.0033 

22.22% 

Table  6.  Turbulence  Model  Performance:  Pressure-derived  Lift  and  Drag  coefficient 
(Ur  =  6m/s) 


Lift  coefficient 

Drag  coefficient 

Experimental 

0.5486 

0.00258 

Predicted 

Error 

Predicted 

Error 

Standard  k-e 

0.5543 

1.03% 

0.0062 

140.31% 

Realisable  k-s 

0.5405 

1.48% 

0.0023 

10.85% 

Standard  k-co 

0.5230 

4.89% 

0.00130 

403.87% 

SST  k-co 

0.5247 

4.55% 

0.0033 

27.91% 

Table  7.  Turbulence  Model  Performance:  Pressure  surface  boundary  layer  parameters  at  93% C 
(Ux  =  3m/s) 


Parameter 

Exp. 

Theo. 

Numerical 

Standard 

k-t 

Realizable 

k-£ 

Standard  k- 

CD 

SST  k-co 

Thickness 

(5) 

23.62 

mm 

30.73 

mm 

78.27 

mm 

27.79mm 

115.82mm 

37.82mm 

Displacement 

thickness 

_ £1 

3.75 

mm 

3.85 

mm 

4.99mm 

3.86mm 

4.94mm 

3.58mm 
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Table  8.  Turbulence  Model  Performance:  Pressure  surface  boundary  layer  parameters  at  93%C 
(Ux  =  6m/ s) 


Parameter 

Exp. 

Theo. 

Numerical 

Standard  k- 

£ 

Realizable 

k~£ 

Standard  k- 

(0 

SSTk-w 

Thickness 

(S) 

24.81 

mm 

26.75m 

m 

92.74mm 

26.44mm 

116.30mm 

38.06mm 

Displacement 

thickness 

(S') 

4.44 

mm 

3.35mm 

4.67mm 

3.74mm 

4.69mm 

3.33mm 

Table  9.  Wall  Treatment  Analysis:  Pressure-derived  Lift  and  Drag  coefficients  (Ux  =  3m/s) 


■  A 

Lift  coefficient 

Drag  coefficient 

Experimental 

0.5520 

0.0027 

Predicted 

Error 

Predicted 

Error 

Standard 

Wall 

Function 

0.4529 

17.95% 

0.0026 

3.70% 

Non¬ 

equilibrium 

Wall 

Function 

0.3933 

28.75% 

0.0036 

33.33% 

Enhanced 

Wall 

Treatment 

0.5302 

4.11% 

0.0024 

11.11% 
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